Bound electron states in the monolayer graphene with the short-range impurities 
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Bound electron states in impure graphene are considered. Short-range perturbations for defect 
and impurities of the types "local chemical potential" and "local gap" are taken into account. 

I. INTRODUCTION 

The Dirac equation is a fundamental base of the relativistic field theory. However, it is an important model in the 
non-relativistic solid state theory as well. Superconductors with d— pairing [T|, the Cohen-Blount two-band model 
of narrow-gap semiconductors [sj, electronic spectrum of the carbon tubes form an incomplete list of the non- 
relativistic applications of this equation. During the last two years extremely much attention was payed to the problem 
of the electronic spectrum of graphene (see for the review ^4]). Two-dimensional structure of it and a presence of the 
cone points in the electronic spectrum make actual a comprehensive study of the external fields effect on the spectrum 
and other characteristics of the electronic states described by the Dirac equation in the 2-1-1 space-time. We consider 
in this work the bound states of the 2-1-1 Dirac equation due to the short-range perturbation. Particular attention 
to this case stems from the effectiveness of short-range scatterers in contrast to the long-range ones: an effect of the 
latter is suppressed by the Klein paradox . Our work takes into account the obvious fact that the Kohn-Luttinger 
matrix elements of the short-range perturbation calculated on the upper and lower band wave functions are not equal 
in a general case. This means that in the perturbed Dirac equation not only the potential but the mass perturbation 
can be present. 

II. PERTURBED DIRAC EQUATION IN (2+l)-SPACE-TIME 

The Dirac equation describing electronic states in graphene reads [3| 

(^-ish cj^d^ - ag (m + 5m) s^^ ?A = (£ - F) V, (1) 

where s is the limiting velocity of the band electrons, cr^ are the Pauli matrices, 2ms^ — Eg is the electronic spectrum 
gap, i}) (r) is the two-component spinor. The spinor structure takes into account the two-band nature. Sm (r) and V{t) 
are the local perturbations of the mass (gap) and the chemical potential. A local mass perturbation can be induced 
by defects in the graphene film or in the substrate [6|. We consider here the delta function model of the perturbation: 

5ra (r) = —bS{r — ro), V{r) = —aS{r — ro), (2) 

where r and rg are respectively the polar coordinate radius and the perturbation radius. Such short-range perturbation 
(and the equivalent form diag{Vi,V2)S{r — rg) with —Vi = ^y^,— V2 = was used in the (3-|-l)-Dirac problem 

for narrow-gap and zero-gap semiconductors in Q. The two-dimensional Dirac problem with the scalar short-range 
perturbation ([2]) (but without the mass perturbation) was considered in \f\ . The obtained there characteristic equation 
for the discrete spectrum energy contains one mistake. We correct it here and take account of the mass perturbation 
Sm (r) . 

Let us present the two-component spinor in the form 

_ expj^iEt) ff^ (r) exp [i (j - 1/2) 0]\ 
U (r) exp [i(i + 1/2)0]]' 
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where j is the pseudospin quantum number; j = ±1/2, ±3/2, . . .. In the opposite to the relativistic theory, this 
quantum number has nothing to do with the real spin and indicates the degeneracy in the biconic Dirac point. The 
upper fj (r) and gj (r) components of the spinor satisfy the equations 

^ + Ig, -{E-m) /, ^{a + b) 6{r - ro)/„ (4) 



- 5 + ^-fj -iE + m) g, = (a - 6) <5(r - ro)ff,. 

These equations have a symmetry; 

fj ^ 9j, -E, j -j 

Let us introduce the function ipj (r) = fj/gi- It satisfies the equation: 

1 



(a + b)Lp] + {a- b) 



Integrating in the vicinity of r = tq 



we obtain the matching condition 



lim 



dipj 



Oj^^^ro-e) {a + b) ip^^ + {a - b) 



+ S{r-ro) =0 



-1, 



arctan (p- 



ja + b 
a — b 



arctan ( ^7 



(5) 



(6) 



(7) 



(8) 
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where ipj = cpj (ro — e) , f'j' = (pj {ro + e) , > b^. The upper and lower component matching condition resulting 
from (l9l) reads 



(10) 



where the matrix A 



cos \/ a?- — b"^ , 

a+b 



'^sinVa^ - 62 



^ sin \/aF~^, cos Va^ h"^ 



(11) 



is orthogonal for 6 = 0. It transmutes into the matrix 



cosh V62 - a2, - V / 1^ sinh V62 - 



'f^sinh \/62 - a2, 



cosh y/b'^ — 



(12) 



when — b"^ < 0. 

The general solution can be found solving the second-order equation obtained by excluding one of the spinor 
components from the equation set ([1]), ([5]) in the domains < r < tq and r > tq : 



h - 0. 



This equation is related to the Bessel one. Its general solution reads 



(13) 



(14) 



3 



where = rri^ — E^, Iv {z) and (z) are the modified Bessel functions. The constant C2 = in the domain 
< r < tq, while Ci = in the domain r > tq. Expressing the t^j-component using ([5]), we can write 



E Ij-1/2 (i^ro) 



V E Ij+1/2 (kvo) ' 



(15) 



E {nr)' 



(16) 



Applying the matching condition © to the expressions (|16p . pS]) we obtain the characteristic equation for the 
bound state energy levels: 



tan (\/^2— 52) 

- &2 



(m - i;) (a - 6) + (a + 6) (m + E) 



Ij+1/2 (uro) Kj+i/2 (uro) 



(17) 



where — b'^ > 0. This equation turns to the characteristic equation obtained in 0], for b = apart from the 
mistakenly omitted terms in the right hand side of ^T7\ . In the opposite case of — 6^ < we have 

Kj-i/2{Kro) Ij-1/2 (kto) 



Kj+i/2{Kro) Ij+1/2 (Kro) 



tanh (V^^ — a^) 
Vb^ - a2 



-(m -£')(&- a) + (5 + a) (rn + £■) 



We write these equations in another form making the symmetry Q manifest: 



1^ [lj-1/2 i^ro) Kj+i/2 {kto) - Kj_i/2 («^o) Ij+1/2 (^''o)] 



(18) 



tan (\/a2 — b'^) 



[(m - i;) (a - 6) Ij+1/2 (uro) Kj+1/2 {kto) + (a + b) (m + E)Ij_i/2 {nro) Kj-1/2 {uro)] , (19) 



K [lj-1/2 (liro) Kj+1/2 (K-ro) - K]-i/2 (k^'o) Ij+1/2 [Kro]] 



tanh - a2) 
- a2 



_-(m - £■) (6 - a) /j+1/2 (k^o) -fsrj+1/2 (^ro) + (6 + a) (m + E)Ij_i/2 {nro) Kj_i/2 (k^q)] (20) 



III. ANALYSIS OF THE CHARACTERISTIC EQUATION AND NUMERICAL RESULTS 

Making use of the Bessel functions limiting forms for small arguments [1] 

I,{z)^{z/2r ^^J^^^ , Ko{z)^-\nz, K^z) ^ (i.) {z/2)-\ 

we can obtain a simple relation describing the asymptotic behaviour of the energy level, where the perturbation power 
approaches zero: 



E = m 



1 - — ^ cxp - 



2r^ V ro (a + b) 



(21) 



4 



where Tc = (in units with h = s = 1), a + b > 0. This resuh conforms the well known general property of the 
two-dimensional quantum systems: a threshold for creation of the bound state is absent; the point a + 6 = is the 
essentially singular point of the function E = E{a + h). One can see that the function E{a) approaches the point 
E = —TO at some large enough value of a > 0. . Making use of the Bessel function asymptotic behaviour Q, 

-1/2 / \ ^^"^ 

Ji, (z) ~ (27rz) ex^pz, Ku{z) ^ y — j exp (-z) , 

and the equation [181 we can see that the function E{h) approaches the point E = Q when ^ is large enough and 

b GO. 

In the Fig. 1 the electron bound state energy is presented as a function of the potential amplitude for the angular 
momentum quantum number j — 1/2, ^ = 1 and 6 = 0. Inspecting this plot one can see that our analytic solution 
()21|1 perfectly approximates approaching of the bound state energy value the upper band bottom, when b approaches 
zero. 

In the Fig. 2 the bound state energy is presented as a function of the mass perturbation amplitude b for a — Q, 
7^ = 1: J = 1/2. 

In the Fig. 3 the electron bound state energy is presented as a function of the potential amplitude for the angular 
momentum quantum number j = 1/2, ^ = 1, and b = —1. We see that the energy dependence on a is non-monotonic 
function, but approaching the upper band bottom takes place similarly to the case of 6 = 0. 



IV. CONCLUSION 



In conclusion, we considered the bound electron states for the two-dimensional Dirac equation with the short- 
range perturbation. The short-range perturbation is approximated by the delta function 5 (r — tq) with different 
amplitudes in the upper and lower bands. We found the characteristic equation for the discrete energy levels. Energy 
levels behaviour in dependence on the perturbation amplitudes was investigated both analytically and numerically. 
The obtained results can be useful for understanding of the graphene electron properties. 
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FIG. 2: Reduced lower electron bound state energy E/m dependence on the mass perturbation amplitude b at a=0. 




FIG. 3: Reduced lower electron bound state energy E/m dependence on the short-range potential amplitude a at b=-l. 



